Quantum spins mixing of spinor Bose-Einstein condensates 
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We examine the internal structure of the ground states 
of a trapped Bose-Einstein condensate in which atoms have 
three internal hyperfine spins. We determine a set of collective 
spin states which minimize the interaction energy between 
condensate atoms. We also examine the internal dynamics 
of an initially spin polarized condensate. The time scale of 
spin-mixing is predicted. 



PACS numbers: 03.75.Fi 



Bose-Einstein condensates (BEC) of atoms with inter- 
nal degrees of freedom are new forms of macroscopically 
coherent matter which exhibit rich quantum structures. 
In the case of BEC with two internal spin states 
theoretical studies have predicted interesting phenomena 
such as quantum entanglement of spins |1 , suppression 
of quantum phase diffusion [Q, interference effects ||. 
Recently, Stamper-Kurn et al. have realized an opti- 
cally trapped BEC in which all three hyperfine states in 
the lowest energy manifold of sodium atoms are involved. 
Such a three-component condensate raises new questions 
regarding the more complex ground state structure [0J|] 
and internal spin dynamics. One of the key features here 
is that there are spin exchange interactions which con- 
stantly mix different condensate spin components while 
the system as a whole remains in the ground state. For 
example two atoms with respective hypefinc spins +1 
and —1 interact and become two atoms with hypefine 
spin 0. Therefore an important problem is to determine 
how atoms organize their spins in the ground state and 
how a spin polarized BEC loses its polarization because 
of spin exchange interactions. 

In this paper we approach the questions using an alge- 
braic method found in quantum optics. We identify the 
fact that the interaction between spin components in a 
BEC is analogous to 4- wave mixing in nonlinear optics. 
However, since the trap is like a matter wave cavity, a 
more appropriate optical analogy is the 4- wave mixing in 
a high finesse cavity (i.e., a cavity QED system). With 
the help of the methods developed in a related cavity 
QED problem [p|JTc|], we are able to study the organiza- 
tion of spins in the condensate ground state. We find 
that there exists a class of quantum superposition states 
which minimize the interaction energy. These quantum 
states are recognized as collective spin states which are 



characterized by strong correlations among different spin 
components, and in some cases we find that the num- 
ber of atoms in individual spin component shows large 
fluctuations. In this paper we also examine the internal 
dynamics of the spin-mixing process arising from the non- 
linear interactions between condensate atoms ]Tl| . For an 
initially spin polarized BEC, we predict the time scale at 
which spins become strongly mixed. 

To begin we consider a dilute gas of trapped bosonic 
atoms with hyperfine spin / = 1. The second quantized 
Hamiltonian of the system is given by (ft, = 1), 

a ^ ' 

/ ¥Jfl^>^ v d 3 x (l) 

where (k — —1,0,1) is the atomic field annihila- 
tion operator associated with atoms in the hyperfine spin 
state |/ = 1, mf — k). The summation indices in (1) run 
through the values —1,0,1. The mass of the atom is 
given by M and the trapping potential Vp is assumed 
to be the same for all three components. The interac- 
tions between atoms are characterized by the coefficients 
QaPfiv which are obtained from the two-body interaction 
model IHH, 

2 F 

U{x l ,x 2 ) = 5(x 1 -x 2 )Y J 9F \F,M F )(F,M F \. 
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(2) 



Here \F, Mp) is the total hyperfine spin state formed by 
two atoms each with spin / = 1, and gp = 4nh 2 ap/M 
with ap is the s— wave scattering length in the F channel. 
The interaction (2) is based on rather general symmetry 
assumptions of the system, because it preserves angular 
momentum and the rotation symmetry in hyperfine spin 
space J8|. The model also makes use of the 5 potential 
approach which has been widely used in one-component 
dilute BEC. 

By expanding the total spin state \F, Mp) in terms 
of basis vectors |/ = 1, m/ = a) ® |/ = 1, mf = we 
obtain the Hamiltonian in the form H = Hs + Ti-A, where 



Hs = J2 d3x ^ 
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Here A s = (170 + 232 )/3 and A a = (32 — .9o)/3 are defined. 
The Hamiltonian TL is written as the sum of a symmet- 
ric part TLs and a non-symmetric part TLa, where TLs 
remains unchanged for any interchange of the spin com- 
ponent indices. 

In this paper we assume that the symmetric interaction 
TLs is strong compared with TLa- This occurs for atoms 
whose scattering lengths in different F channels have ap- 
proximately same values such that |A S | ^> |A a |. Recent 
estimations have indicated that sodium and rubidium 
atoms indeed have such a property. With the symmet- 
ric TLs being the dominant Hamiltonian, the condensate 
wavefunctions 4> K (x) (k = 0, ±1) for each spin compo- 
nent are approximately described by the same wavefunc- 
tion <fi(x), i.e., <f> K (x) = 4>(x), which is defined by the 
Gross-Pitaevskii Equation through TLs , 



.Yl + V T + Xs N\q>\ 2 



(5) 



where e is the mean field energy or the chemical potential. 

Under the condition that atoms in different spin states 
are described by the same wavefunction, we can approx- 
imate field operators at the zero temperature by, 



K = 0,±1 



(0) 



Here a K is the annihilation operator associated with the 
condensate mode, and it satisfies the usual commutation 
relation [d K ,d 7 ] = and = § K1 , Using (5) and 

(6), TLs and TLa have leading parts H s and H a respec- 
tively, 

TLs~£N- X' s N(N -1) = H s (7) 
TLa ~ X' a (ji\a\aiai + alja^a-xfl-i + 2djdJdido 
+ 2a! 1 a ( ! ) d_iao — 2d|dL 1 did-i + 2djd ( ! ) did_i 
+ 2d|dl 1 d d ) = H A (8) 



Here 2A- = Aj J |</>(x)| 4 d 3 x (i — s,a), and N = a\a\ + 

d^do + SL1S-1 is the total number of atoms in the con- 
densate. 

Our goal is to find the quantum states that mini- 
mize the energy H s + H a . Since H s is a function of 
N only, H s is a constant operator for a fixed number 



of atoms. Therefore it is sufficient to look for the ground 
state of H a . It is quite remarkable that a similar struc- 
ture of H a also appeared in nonlinear wave- mix ing pro- 
cesses in cavity QED [||. We follow Refs. HJnJ and 
identify the algebraic structure of the system. We no- 



tice that the operators L_ = \[2 
\/2 ^dpfii + d^do^ and L z = 



a\a 



L + = 



L d_i — djdi ) obey an- 



gular momentum commutation relations: 



2L, and 



L z ,La 



= ±L±. In other words, the operators 



L + , L_ can be interpreted as raising and lowering oper- 
ators of a kind of 'orbital angular momentum', and L z is 
the 'z-component' in the standard notations. From the 
theory of angular momentum, L 2 and L z have a complete 
set of common eigenvectors \l,mi) defined by 



L 2 \l,mi) =l(l + l)\l,mi) 
L z \l,mi) = mi \l,mi) 



(9) 
(10) 



where mi — 0,±1,±2, ...,±Z. For a given total number 
of atoms N, the allowable values of I are / = 0, 2, 4, ...N 
if N is even, and I = 1,3, 5, N if N is odd. 

With the help of the angular momentum operators, H a 
takes a very simple form, 



H a = K (t 2 - 2N) 



(11) 



This is the main result of the paper because the energy 
spectrum of H a is now solved. Eq. (11) indicates that 
\l,mi) are eigenstates of H a with the energy Ef 



E? = X' a [1(1 + 1) - 27V] 



(12) 



The lowest energy state of H a depends on the sign of A„ . 
In the following we discuss two cases: (I) \' a > and (II) 

K < 0- 

(I) X' a > 0: In this case \l = 0,m/ = 0) is the ground 
state of H a - Using the Fock states \m,no,ri-i) defined 
by the number operators hj = djdj for the three spin 
components (i.e., hj \n\, no, n_i) = Uj \ni, no, ri-i)), 
\l — 0, mi =0) has the form 

[N/2] 

\l = 0,mi = 0) = A k \k,N- 2k, k) (13) 

fc=0 

where the amplitudes obey the recursion relation 



N-2k + 2. 
w N -2k + l 



(14) 



We see that the state \l = 0,m; = 0) is a quantum su- 
perposition of a chain of Fock states \k, N— 2k, k) in 
which the numbers of atoms in the spins 1 and —1 are 
equal. We stress that such a quantum state is a collective 
spin state which cannot be expressed as product states 
of individual atoms. The amplitudes Ak are arranged 



2 



in such a way that the interaction energy H a is almost 
completely cancelled. This can be seen from the disap- 
pearance of ./V 2 dependence in the energy of H a . It is 
not difficult to show that for the state (13), the aver- 
age numbers of atoms in each component are all equal, 
i.e., (no) = (hi) = (h-i) = N/3. Since are almost 
uniformly distributed (see Fig. la), there are large fluc- 
tuations of particle numbers in individual components 
although the total particle number N is fixed. More pre- 
cisely, we find that (An ) « 2N/y/E for N » 1, i.e., a 
super-Poisson distribution. Our further calculations indi- 
cate that super-Poisson distribution of particle numbers 
are a common feature for low energy eigenstates of H a 
when > 0. 

(II) X' a < 0: In this case H a has 2N+1 degenerate ground 
states given by |Z = N, mi) where mi — 0, ±1, ±2, ±N. 
The energy (12) of these states is X' a N(N — 1), and the 
general form of \l — N, mi) is given by 

\l = N,mt) =Y. B( k n ' ) \k,N-2k-mi,k + mi). (15) 
k 

Here the summation index k runs over all physical Fock 
states \k, N ~ 2k — mi, k + mi) (i.e., those with non- 
negative numbers in each component). The simplest case 
of (15) is \l = N,mi = -N) = |JV, 0, 0>, and with this we 
can construct the amplitudes B^ 1 by repeatedly ap- 
plying the raising operator L + . To give an illustration, 
we plot in Fig. lb the Fock state amplitudes -B^" 1 '' for 
several mi's. We see that B ( ™ l) has a narrow distribu- 
tion which indicates well defined particle numbers in each 
spin components. It is interesting that all the degenerate 
states (15) have sub-Poisson number fluctuations in each 
spin component (see the inset of Fig. lb). This feature 
is just the opposite to the previous case A' a > 0. 

Finally, let us look at the spin-mixing dynamics of an 
initially spin-polarized condensate in which all atoms in 
the condensate are prepared in the spin state at t = 0, 
i.e., \ip(0)) = |0, N, 0). In this case two atoms in the spin 
state can be converted into one atom in the spin 1 state 
and the other in the spin —1 state. Using the H s + H a 
as our approximate Hamiltonian, the system at time t is 
given by, 

N 

= e- ie " (t) Cie~ tX '« l{l+1)t \l,mi = 0) (16) 

1=0 

where d = (l,rni = 0\ 0,N,Q) and 9 N {t) = (sN + 
X' S N(N — l))t. In Fig. 2 we present the time depen- 
dence of the particle number in the spin-0 component for 
N = 10 2 , 10 3 , 10 4 cases. We see that the number of atoms 
in the spin-0 component becomes steady at (no) = N/2 
after a time t c , 



This is the time scale for the spin-mixing process purely 
due to the nonlinear interaction between condensate 
atoms [Q. In the Thomas- Fermi (large N) limit, we 
find that for a spherial harmonic trap, (g 2 — 9o)tc ~ 

5-liV 1 / 10 ((g + 2g 2 )/Mu) 2 ) 3 ^ where u is trap frequency. 
Therefore t c becomes quite insensitive to N in the 
Thomas-Fermi limit. To give a realistic example, for a 
sodium condensate with N — 10 4 and uu = 27rx370Hz, 
we find that t c is about 0.4 seconds, assuming a 2 ~ 2.6nm 
and ao w 2.3nm. 

We remark that the spin mixing dynamics can be quite 
different for different initial conditions. In Fig. 3, we 
give an example for the case when all three components 
initially have the same atom numbers, i.e., \ip(0)) — 
| AT/3, N/3, N/3). It is quite surprising that the parti- 
cle number executes fast oscillations with a frequency 
of the order of X' a N, and then the system suddenly be- 
comes steady. This interesting behavior indicates that 
there are complex quantum dynamics governed by the 
nonlinear interaction H a . In fact since H a has a discrete 
spectrum, quantum recurrence or revival is expected in 
a much longer time scale (which is typically of the or- 
der of n/X' a ). It is worth further exploring the quantum 
dynamics in the context of either BEC or cavity QED. 

To conclude, we have examined the spin-mixing inter- 
action of a Bose-Einstein condensate with three internal 
spin components. It is quite remarkable that the model 
interaction (2) (which is based on general symmetry as- 
sumptions) can lead to a simple algebric representation, 
and from which we can construct the collective spin states 
which minimize the interaction energy among condensate 
atoms. These collective states exhibit spin correlations 
and characteristic particle number fluctuations which de- 
pend crucially on the sign of X' a . We have also investi- 
gated the spin mixing dynamics due to the nonlinear in- 
teraction between condensate atoms. The time scale of 
mixing for an initially spin polarized system is identified. 
This study provides a theoretical treatment of the struc- 
ture and dynamics of spinor BEC. However, our analysis 
are limited to interaction between condensate atoms, it 
remains to be answered that how non-condensate atoms 
will decohere the condensate structure. Further work 
along this direction would be necessary. 
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FIG. 2. Time dependence of average number of atoms in 
the spin state normalized by the total number of atoms TV. 
The initial state of the system is |^>(0)} = 0,7V, 0). We show 
three cases with TV = 10 2 , 10 3 , 10 4 . 
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FIG. 1. Amplitudes of Fock states associated with the 
ground states of H a for N — 10 3 atoms: (a) > 0, (b) 
X' a < 0. The inset in (b) shows the normalized number fluc- 
tuations Aj = (Arij) 2 I (hj) in the three spin components 
(j — 0, ±1) as a function of mi. Sub-Poisson distributions are 
defined by Aj < 1. 



FIG. 3. Time dependence of average number of atoms 
in the spin state. The initial state of the system is 
|V(0)> = \N/3,N/3,N/3), where N = 300. 
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